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Following the original approach of Maxwell-Boltzniann(MB), we derive a 4- velocity distribution function for 
the relativistic ideal gas. This distribution function perfectly reduces to original MB distribution in the non- 
relativistic limit. We express the relativistic equation of state(EOS), p — po = (7— in the two equations: 

P — PofWi ^iid p — pog(X), where A is a parameter related to the kinetic energy, hence the temperature, of 
the gas. In the both extreme limits, they give correct EOS: p ~ 3p in the ultra-relativistic, and p — po — 
in the non- relativistic regime. Using these equations the adiabatic index 7 (=^) and the sound speed as are 
calculated as a function of A. They also satisfy the inequalities: | < 7 < | and < perfectly. 
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I. INTRODUCTION 

The thermodynamics of matter plays an important 
roll in the study of its fluid dynamics. The solution 
of hydrodynamic equations, both in the relativistic and 
non-relativistic formalism, require a knowledge about the 
equation of state (EOS) of the system. An EOS is essen- 
tially a relation between the macroscopic quantities e.g. 
the pressure p, the total mass-energy density p, and rest 
mass energy density po- One could arrive into such a re- 
lation from the knowledge of the probability distribution 
of various microscopic quantities (e.g molecular velocity, 
energy). For the relativistic gas, such a study has been 
done using the formalism of statistical mechanics in a 
Lorentz covariant framework (See [l| and ref. therein, 

II, II, 11) • Distribution function of co-ordinate velocity 
(^) and the EOS can be derived using this formalism. 
Several authors tried to use relativistic EOS to solve the 
relativistic hydrodynamics (RHD) equation in various as- 
trophysical problems (5]). But a direct application (Q) 
of the exact relativistic EOS turns out to be a less feasible 
method for numerical computing because of its complex- 
ity. Thus in the popular studies of numerical RHD in 
astrophysical problems (see review 7]), people prefer to 
use various alternative models of EOS which can repro- 
duce various features of exact EOS with a good approxi- 
mation 0). All these models ([l3l) are proposed in 
an ad-hoc manner and are not derived consistently from 
the first principle. 

The present status of RHD study motivates us to ana- 
lyze this problem in a model independent way in partic- 
ular, we wish to find a simpler form of EOS which, apart 
from being computer friendly, must also be derived from 



very first principle of kinetic theory. At this point, it 
would be interesting to see whether one could derive a 
probability distribution function of 4-velocity following 
the basic methods of probability theory and assumptions 
of isotropy as originally had been done by MB in the case 
of non-relativistic gas. This approach, to the best of our 
knowledge, is yet unexplored and one could examine it 
atleast for the sake of conceptual completeness. An EOS 
could be derived then using this distribution function. In 
this article, we follow the original line of arguments of MB 
to find a distribution function of 4-velocity (u* = ^), 

instead of co-ordinate velocity ^ , for a relativistic ideal 
gas. We then use these Maxwell like 4-velocity distribu- 
tion function to derive a form of EOS. 

This paper is organized as follows. In §11, we derive 
Maxwell like 4-velocity distribution function for relativis- 
tic gas. In §111, we apply it to find the EOS for rela- 
tivistic gas. In §IV, we investigate the behavior of some 
physical quantities (e.g. 7 and sound speed etc.). The 
extreme limits: ultra-relativistic and non-relativistic are 
discussed. Finally concluding remarks are drawn in §V. 



II. FOUR VELOCITY DISTRIBUTION OF A 
RELATIVISTIC IDEAL GAS 



We consider a gas consisting of one species of non- 
interacting particles. The temperature of the gas is so 
high that the average kinetic energy of a constituent par- 
ticle is comparable to its rest mass energy. In other 
words, the thermal energy of the gas is comparable with 
its rest mass energy. The co-ordinate velocity and the 
4-velocity of a gas particle are defined respectively as 
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dt 



and 



dr 



where rfr^ ~ c?dt^ — {dxY is the 



length element in the Minkowski space-time, u'^'s are 
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related to the co-ordinate velocity as 



(1) 



where v'^ = (!,«)• The molecules of the gas interacts 
only through elastic collision so that all of their energy 
is kinetic energy. In a Lorentz frame, where the center 
of mass of the container is at rest, the distribution of 
molecular speed is isotropic. The choice of such a frame 
is unique upto a three dimensional rotation. In different 
frames in which the container is at rest, takes same 
value and u*s are transformed by a 3-D Euclidean rota- 
tion i. e they transform like a components of an ordinary 
three dimensional vector in non-relativistic case. Accord- 
ing to special relativity 



(2) 



where — is the co-ordinate velocity. In the non- 
relativistic case (where Vx,Vy, Vz are not constraint by 
the above inequality) , we write the probability of finding 
the particle simultaneously in the velocity range to 
Vx + dvx, Vy to Vy + dvy, aud, Vz to Vz + dvz as the 
product of the individual probabilities: 

p{vx,Vy,Vz)dvxdvydvz ^ f{vx)f(.Vy)f{vz)dvxdvydvz (3) 

where f{vi)dvi represents the probability of finding the 
particle in the velocity component range Vi to Vi + dvi . 
The functional form of '/' is same for all the components 
due to isotropy. This assumptions ^ is no longer valid 
in relativistic case as the co-ordinate velocity component 
are no longer independent due the above inequality and 
can only vary from — c to c. This problem can be cir- 
cumvented if we seek a formula for the probability distri- 
bution of 4-velocity instead of coordinate velocity. The 
4-velocity components are related through the equation 



(4) 



But, as u — ^ c, u*^ — > oo and therefore, each of the 
three spacial components (u\ u^, u^) of the 4-velocity can 
take values from — cxd to cx) independently, irrespective 
of the values of other two components. Therefore, just 
as the 3- velocity components in the non-relativistic case, 
the three spacial components of 4- velocities have no con- 
straint among themselves. In the chosen Lorentz frame, 
the center of mass of the gas system is at rest. The av- 
erage 4-momentum of the total gas system is zero in all 
directions. This implies that the probability distribution 
function of the three spacial components of the 4-velocity 
is isotropic just as the probability distribution function 
for velocities in non-relativistic case. Let the probability 
of finding a gas molecule in between to + du^ be 



p,{u')du' = F,{u')du\ {i ^ 1, 2, 3). 



(5) 



The assumption of isotropy demands that the functional 
from of Fi is same for all values of i, i.e. Fi = F2 = -F3. 



As already explained above the three spatial components 
of 4-velocity m^, m^, are not constrained, hence the 
probability that a gas particle is found simultaneously 
within the 4-velocity range to +du^ , to u^ + du^, 
and to + du^ is 

p{u\u^,u^)du^du^du^ = F{u^)F{u^)F{u^)du^du^du^ . (6) 

But, on account of isotropy, p{u^, w^, u^) is invariant un- 
der a three dimensional rotation and must be a function 

of (7/1)2 ^ ^^2-)2 ^ ^^3)2 ^ ^y)2^ rpj^jg gjygg 



F{u^)F{u^)F{u'^) ^ iPiiu)^). 



(7) 



Differentiating EqjT] partially, w.r.t u^, , and v? , and 
dividing \yYF{u^)F{v?)F{u^), one arrives at the equation 



F'{u^) _ F'ju"^) _ F'ju^) _ i^'ju) 
uiF(ui) ~ v?F{v?) ~ u^F{u^) ~ uiP 



= -2A (8) 



where A is a constant which is positive because of the 
convergence requirement of F. The Eq.® is solved to 

give: F{u^) — Ae^^^'^ ^ and therefore. 



F(ui)F(ii2)F(u3) = A^e-^t")' 



(9) 



Integrating this probability over spherical shell of radius 
u and u + du, one finds the total probability of finding 
the particle between the speed range u and u + du, as 



F{u)du = AnA^u^e-^'' du. (10) 

The normalization condition: F{u)du = 1 relates A 

with A as A = \ -. The constant A is related to the 
average energy of the gas particle as 



u^F{u)du ^ AttA^ / (1 + — ) e 



-Xu' 



du. 



As in the case of MB, the parameter A in this case, is 
also related to the average kinetic energy {<Ek>) of the 
gas molecule, however not in a simple manner. If m is 
the rest mass of the gas molecule then 



<Ek{X)>^ mc2[<7iO> -1] 



-Xu' 



du 



mc 



(11) 



The above equation gives the physical meaning of the 
constant A as it relates this constant to the physically 
measurable quantity < Ek>, i.e. A is a measure of the 
temperature T of an ideal gas. < Ek >, and T are 
monotonically decreasing function of A. Therefore, the 
non-relativistic limit is achieved in large A limit. The 
function F{u) has non-negligible value only in the region 
where Xu^ ^ 1 or, u^ ^ j. For large A, this implies 

= v'^/{l - ^) < 1, i.e. ^ < 1. This gives « v\ 

therefore, F{u)du = AnA^v'^e'^^ dv, is the velocity dis- 
tribution function in the non-relativistic case. 
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III. THE RELATIVISTIC EOS 

In the general case of a relativistic gas, one can relate 
p, p, po and the adiabatic index 7 (= — ) as 



P= (7- l)(p-Po)- 



(12) 



It is easy to show that for a cool non-relativistic gas 7 
becom es | , and for a extreme relativistic gas of photon 
7=1 [llj . Therefore, the value of 7 in general depends 
on the relativistic nature of the gas, i.e. the dominance of 
the kinetic energy density over the rest mass density pq. 
To calculate the dependence of 7 on p, we first expressed 
p, and p, in terms of po- This can be done by considering 
that the pressure p is the momentum flux through unit 
area averaged over all the molecules of the fluid-blob in 
the co-moving frame. We take a surface element of the 
fluid blob oriented to any one axis (z*'' axis say,) the 
pressure p is p — n < p^v^ >. Here n is the number 
density of the molecule, is the i*'' component of the 
4-momentum p'^ = mu'^. Using the relation = and 



(" ) 



> = < 



taking < ^"o >=< — ^-^ — >^ — 

o 

isotropy, we write the pressure as 



liny 

3 uO 



> due to the 



P^\po{<y°> - <^>) = Po9{\)- (13) 

Here and henceforth we work with c — 1. The total 
energy density p is related to rest mass density po ^-s 

p^nm <u*'(A)>= Po <u"(A)>= po/(A), (14) 

where /(A) and g{X) are related to the modified Bessel 
functions of second kind A'„(A) through the 

<7/°(A)>=47rA3/yTT^u2g-Au^^^^ (^5) 
Jo A 2 



<—^>=^TrA-^ 



~mO(A)' 



VT 



-^o(^)].(16) 



The properties of -ftr„(A) functions and its derivatives 
(require later) are well known in mathematics ([l2])- The 
value of AttA^ is given before using the normalization 
condition. 
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FIG. 1. The variation of specific internal energy {e) and spe- 
cific enthalpy (h) are shown as a function of A. 



as function of A only. To find the expression for the 
sound speed, one needs to start from the definition aj. = 
(fp)sj together with the entropy equation Tds = and 
P = p{p,Po)- Computing the quantities (^)s, (|i)po> 



il 



-)a, and simplifying one finally yields 
1, dp 



'^dp'"'' 



1- 



op . 



+ f(^)A-(l + !7)(/ 



h dpo' 



■), (18) 



in terms of A completely. Here the prime ('/') denotes 
derivatives w.r.t A. To check the overall variation of 7 and 
tts w.r.t A, we plot these functions in Figl5]and in FiglJl 
In Figl2l we see that there is no applicable variation of 7 
throughout. However there is a sharp decay of 7 values 
very close to zero. This indicates that EOS is mostly non- 
relativistic in nature and becomes relativistic to ultra- 
relativistic at the end phase of evolution. On the other 
hand, the deviation of the sound speed in Fig|31 from 
non-relativistic value is prominent even in the relatively 
higher A values. 



IV. EXTREME LIMITS: ULTRA-RELATIVISTIC AND 
NON-RELATIVISTIC 



A. The Specific Heat Ratio and the Sound Speed 

Using the EOS ([131 [M]), one could to find the varia- 
tion of 7 w.r.t A for a relativistic gas in many astrophys- 
ical applications. Therefore, our next job is to express 
the thermodynamic variables in terms of A alone. Using 
thermodynamic definitions, we obtain the specific inter- 
nal energy e ~ ^^f^ = / ~ l^ ^-nd the specific enthalpy 
h = = / + 5 of the gas. Variation of e, h are plotted 
in Fig[T] as a function of A. Further, simultaneous solu- 
tions of three equations: p = (7 — l)(p — Po), P = Po/(^)i 
and p = pQg{X), provide the specific heat ratio 



7=(/-l)-i(/i-l), 



(17) 



It is important to show that this expression for p gives 
correct results for non-relativistic and ultra-relativistic 
limit. In the non-relativistic limit v << c. Therefore, 
■° ^ (1 + 5f^) and ^ « (1 - ^^). Therefore, 



P 



u 

■ipo >. Now, the average kinetic energy density is 



po < Ek >= po(< w° > -1 



ipo < >. 



Therefore, 



p — Pq = ^po <v >. Using above results, we find 



PO = 



(19) 



an exact EOS in the non-relativistic limit in which 7 = §■ 
For an ultra-relativistic gas of photon u — >• c and m — 0. 

Ephoton remains finite. Therefore, for 
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FIG. 2. A plot of the sound speed as vs. A is shown in figure. 
The inequality, Os < l/y/S satisfies when A — >■ 0. 
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FIG. 3. The adiabatic index 7 (| < 7 < |) is plotted w.r.t A. 



a photon gas u° — > oo and — > 0. This gives po = 0, 
p = nEphoton and p = \nEphoton- Therefore, in the 
ultra-relativistic limit, EOS and 7 respectively becomes 



P=\p, 



(20) 



For a gas of massive particle, the values of the 7 and 
fls in the extreme relativistic regime can be obtain from 
Eq.dlTl) and Eq.dlHl) in the limit A ^ 0. The integral 
in Eq. ([T5| & ([T6l) and their derivatives w.r.t A, take the 
value for A — s- (by putting ^v? = y) 



2A- 



-,<u''> '=- 



A- 



, <- 



1 



2A5 



1 



:> =■ 



A- 



an exact match, whereas, in the non-relativistic regime, 
EOS correctly reduces to p — po = |p, which implies 
7 = |. The results we obtain are also in well agreement 
with the known results. Our theoretical investigation 
indicates that EOS evolves abruptly at the end phase of 
ultra-relativistic regime. Following the prescription of 
new EOS, it is easy to express the thermodynamic and 
hydrodynamic variables in terms of A alone. In the next 
stage, one just need to estimate the variation of A w.r.t 
the radial distance as to find an interesting solutions in 
different physical systems in interest. We have studied 
one such solutions in the accretion around black hole 
which we wish communicated soon. 



Substituting above values in the Eq. p7|) . & (|18p . we find 



lim 7 = - 



and lim a? = — . 

A^o 3 



(21) 



Similarly, it can be shown that in the large A limit, 
7 = |. We emphasis these points in Figl2l where y-axis 



are set at y — ^ and |. We see, 7 
holds in our case. Also, in the Fig[ 
speed as < satisfies. 



< 



7 < I perfectly 
we find the sound 



V. CONCLUSION 

In this communication we express the relativistic 
EOS in a new form through two parametric equations: 
p = pof{X), and p = pog{X), where A is a parameter 
related to kinetic energy of the gas. The EOS is 
obtained by using a 4-velocity distribution function 
which we derived by applying the basic principles of 
probability theory with the assumptions of isotropy. 
In the ultra-relativistic regime, these new equations 
perfectly produces well-known results, namely, the EOS: 
p — 3p, the value of 7 = | , and sound speed < ^ have 
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